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Abstract 

We argue that to solve the foundational problems of quantum theory one has to first understand 
what it means to quantize a classical system. We then propose a quantization method based on 
replacement of deterministic c-numbers by stochastically-parameterized c-numbers. Unlike canon- 
ical quantization, the method is free from operator ordering ambiguity and the resulting quantum 
system has a straightforward interpretation as statistical modification of ensemble of classical tra- 
jectories. We then develop measurement without wave function collapse a la pilot-wave theory and 
point out new testable predictions. 
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I. INTRODUCTION 



While pragmatically quantum theory has shown spectacular successes, its foundation 
with respect to quantum-classical correspondence and measurement problem, despite many 
attempts, still resist unambiguous explanation. We believe that these difficulties originate 
from the (physical) ambiguity of the quantization scheme through which many successful 
quantum systems are obtained from the corresponding classical systems. In other words, 
to solve the foundational problems of quantum theory, we have to first understand what is 
meant by quantizing a classical system. In the canonical quantization, one first writes the 
classical dynamical equation in Cartesian coordinate. Quantization is then done by promot- 
ing the pair of canonical conjugates variables into the corresponding Hermitian operators, 
o i — y 6, and Poisson bracket is replaced by the commutator: {°,*} i— > [o,i]/ih, where h is 
the reduced Planck constant. This procedure is usually said as replacement of commuting 
c-number (classical number) by non-commuting q-number (quantum number or Hermitian 
operators) . 

Such direct substitution rule however implies formal, conceptual, and foundational prob- 
lems. First, given a classical quantity, the above rule in general leads to infinitely many 
different alternatives of Hermitian operators due to operator ordering ambiguity. Second, 
the physical meaning of the resulting quantum systems can not be deduced "directly" from 
the quantization processes: the quantization does not tell us the physical meaning of the 
Schrodinger equation, nor it offers an explanation to the physical origin of Planck constant. 
In other words, the quantum-classical correspondence is not physically transparent. Hence, 
unlike the original classical system, we need to further introduce a physical interpretation 
to the resulting quantum system. Unfortunately, there are many physical interpretations 
with the same empirical prediction: standard Copenhagen interpretation, pilot-wave theory, 
many worlds etc. There is then the foundational problem that that some of the existing 
interpretations suffer from the infamous measurement problem In the present paper, 
we shall discuss a new method of quantization which are free from the above mentioned 
problems. 
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II. QUANTUM FLUCTUATIONS AS STATISTICAL MODIFICATION OF CLAS- 
SICAL ENSEMBLE PARAMETERIZED BY UNBIASED HIDDEN RANDOM 
VARIABLE 



A. General formalism 

Let us consider a classical dynamics with N degree of freedom whose configuration co- 
ordinate is q — (qi, . . . , qw). Generalization to infinite degree of freedom is straightforward 
at least formally. In classical mechanics, all the dynamical information is contained in the 
Lagrangian function L = L(q,q), where q = dq/dt, is the velocity. For simplicity, let us 
assume that the Lagrangian is non-singular. From the Hamilton principle, one then has the 
Euler-Lagrange equation (d/dt)(dL/dqi) — dL/dqi = 0, i = 1, . . . , N. To go to the Hamil- 
tonian formalism, first, the canonical momentum is defined as p. = dL/dqi, i = 1, . . . , N, 
which, due to non- singularity of the Lagrangian, can be inverted to give <ji = <ji(q,p]t), 
i = 1, . . . ,N. The classical Hamiltonian is defined as H_{<liP) = P^i^iP'-, t) — L{q,P', t). The 
Euler-Lagrange equation can then be put into the Hamilton equation 

* = ^, v_=-^ « = i,...,tf. (1) 

Now let us consider a congruence of curves, {q(t; u) , p(t; v)} , satisfying the Hamilton 
equation, where u = {ui, . . . , u^} and v = {vi, . . . , vjy} are parameters of the congruences. 
Namely, each pair of the value of (u, v) corresponds to a single member of the congruence. 
Let us then assume that there is a one-parameter family of hypersurfaces S_(q; t) = r, where 
S_ is a real-valued function and r is a parameter, satisfying the following relation: 

Pi = d qi S. (2) 

Then, it can be shown that S(q; t) also satisfies the following Hamilton- Jacobi equation j^J: 

d t S + H(q,d q S) = 0. (3) 

Conversely, defining p as in Eq. (T5]), the Hamilton- Jacobi equation of ((3]) can be shown to 
lead to the Hamilton equation of (OQ). 

Now let us consider an ensemble of copies of the system. Then, if the probability density 
of the position of the ensemble of trajectories in configuration space is denoted by p{q',t), 
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then it must also satisfy the following continuity equation: 

d t p + d q -(q(S)p) = 0, (4) 

where the functional form of the velocity q with respect to S_ is determined by substituting 
Eq. ([2]) into the left equation of ([I]) 

m) = 1 = , i = i,...,n. (5) 

Op. p=dgS 

Hence, the dynamics and statistics of the ensemble of trajectories is given by solving Eqs. 
©, ® and ©. 

Let us develop a general scheme to modify the above classical dynamics of the ensemble of 
trajectories To do this, let us introduce two real- valued functions A; t) and S(q, A; t) 
where A in non-vanishing hidden random variable, A ^ 0. A; t) is the joint-probability 
density of the fluctuations of q and A so that the marginal probability densities are given by 

p(q;t) = J d\Q, P(A) = J dqQ, (6) 

where we have assumed that the probability density of A is stationary. 

Now let us postulate the following rule of replacement of (deterministic) c-number by 
(stochastic) c-number to be applied to Eqs. ([3]) and (jlj) governing the dynamics of the 
ensemble of trajectories: 

p^Q, 9g S^d q S + ^, d t S^d t S+^ + ±d g .q(S), (7) 

where the functional form of q{S) in the third line is determined by the classical Hamiltonian 
as in Eq. (^j. That is we just substitute S in Eq. ([5]) with S. Let us show that the above rule 
of replacement has a consistent classical correspondence. First, expanding as AF = dfFAt+ 
d q F-Aq, using the last two equations of © one gets AS ^ AS+(\/2)(Att/tt + d q -q(S)At). 
In the limit of S — > S_, the second term on right hand side has to be vanishing to give 
AQ + Qd q ■ q(S_)At = 0. This is just the classical continuity equation so that one has 
p = J dXfl — > p. 

The next question is then what is the statistics of A. We shall show in the next subsection 
by taking a concrete example that the above modification of classical dynamics of ensemble 
of trajectories leads to the Schrodinger equation with unique quantum Hamiltonian if the 
probability density of A is given by 

P {X ) = 1 -5(\-h) + ±5(\ + h). (8) 



Namely A is a binary unbiased random variable which can take values ±h. Extension to a 
more general continuous hidden random variable satisfying the unbiased condition 

P(X) = P(-\), (9) 

will be given in Section 4, suggesting a very small yet finite corrections to the prediction of 
quantum mechanics. 



B. Particle in external potentials 

Let us apply the above general formalism to a single particle subjected to external po- 
tentials. The Lagrangian is given by L = ^gij(q)q l q j + A ■ q — V(q), where A = (Ai, A 2 , A 3 ) 
and V are the vector and scalar potentials respectively, gij(q) is an invertible matrix which 
might depend on q (for example in the case of particle with position-dependent mass widely 
used in solid state physics), and summation over repeated indices is implied. Writing the 
inverse of as g^ so that gikg kl = the momentum is related to the velocity as 

c? = g^{ P _ r A ). (10) 

so that the classical Hamiltonian takes the following form: 

K(q,p) = - A^p. - A,) + V. (11) 

The Hamilton- Jacobi equation of ([3]) then reads 

d t S + ^-{d qi S - Ai){d qj S - Aj) + V = 0. (12) 

On the other hand, putting Eq. d2J) into Eq. (fTU]) . the functional relation between q and S_ 
is given by 

= 9 ij (d q .S — Aj). (13) 
Inserting this into Eq. fl3j), one thus has 

d.p + d^id.-S-A^p) =0. (14) 

Hence, the dynamics and statistics of the classical trajectories are determined by solving 
Eqs. (H2D, (H3H and (HU). On the other hand, using Eq. (H3J, Eq. © becomes 



p^ Qj dq s^d q S + ^, d t S^ d t S + ^ + ^d qi9 v(d q] S - A 3 ), (15) 



Now let us apply the rule of replacement of Eq. (fl5l) to Eqs. ( TT2|) and ffT4"|) . First, 
imposing the first two equations of Eq. ffl5l) into Eq. ffT4l) . one has 



d t n + ^ (g^s- Aj)n) + ^d qi ( g »d q .n) = o. (ie) 

On the other hand, inserting the last two equations of ffl5|) into Eq. ffl2l) . one has, after 
arrangement, 

^ + - A t) (^5 - ^ + y - £ (^'^ + 

+^(ftn + a ft (^"(^-^)n) + ^a ft 07%n)) =o, (17) 

where we have defined R = Vu and used the following identity 

i a g . ^ = i d qi d qj n a qi d qj r 

4 n 2 2 n r ' 1 J 

Inserting Eq. ( Fl6l) into Eq. ( flTl) . the second line is vanishing to give 

^ + 9 ^{d qi S - A t )(d qj S - Aj) +V- £ (g- 9 -^- + d qi g^) = 0. (19) 

We have thus pair of coupled equations ffl6|) and ffl9|) which are parameterized by A. 
Now let us assume that A; t) has the following symmetry: 

n(q,\;t) = n(q,-\;t), (20) 

so that the probability density of A is unbiased -P(A) = J dqQ(q, A; t) = P(— A). In this case, 
S(q, A; t) and 5*(g, —A; t) satisfy the same differential equation of (fT9~l) : namely the last term of 
Eq. f|T9l is not sensitive to the signs of A. Hence, if initially one has S(q, A; 0) = S(q, —A; 0), 
then the symmetry will be preserved for all the time 

S(q,X;t) = S(q,-X;t). (21) 

These properties can then be used to eliminate the last term of Eq. (fl6|) : taking the case 
when A is positive add to it the case when A is negative and divided by two, one gets 

d t Q + d qi (g ij (d qj S - Ajri) = 0. (22) 

We have thus pair of coupled equations f|T9|) and fT22|) which are still parameterized by A. 
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Next, since A is assumed to be non- vanishing, one can define the following complex- valued 
function: 

y(q,\;t) = Rexp^Sy (23) 

It differs from the Madelung transformation in which S is divided by |A| instead of h. 
Equations (fT9j) and (122]) can thus be recast into the following modified Schrodinger equation: 

i\\\dt* = ^(-i\X\d qi - ^)^'( g )(-i|A|9 3 . - Aj)V + V^. (24) 

Here we have assumed that the fluctuations of A in space and time are ignorable as compared 
to that of S. Let us proceed to assume that Q is separable Q(q,\;t) = p(q;t)P(\), where 
P(A) takes the form given by Eq. (jSj). In this case, Eq. (J2"ll) reduces into the celebrated 
Schrodinger equation 

ihd t ^ Q {q; t) = HV Q (q; t), with ^ Q (q; t) = Jpe* 8 **®, and S Q (q; t) = S(q, ±h; t), (25) 

where the quantum Hamiltonian H is given by 

H = hip* - A^gv^ipj - Aj) + V, with p t = -m qi . (26) 

From Eq. ff25l) we know that the Born's statistical interpretation of wave function is valid 
by construction for all time, p{q;t) = \^Q(q; t)\ 2 . 

Let us mention that there are many approaches to derive Schrodinger equation with 
quantum Hamiltonian of the type of Eq. (1261) |4|. The advantage of our approach is that 
it can be applied directly as soon as the classical Hamiltonian (or Lagrangian) is given, if 
a solution exists. Our method of deriving the Schrodinger equation can thus be regarded 
as to provide a method of quantization of a classical Hamiltonian. For the case gij = mdij 
where m is constant so that g l i = 5ij/m, we regain the result of canonical quantization. 
For the case where g^ depends on q, then in contrast to canonical quantization which leads 
to infinite different alternatives with different ordering of operators, our model imposes a 
unique ordering. The same ordering is also obtained in Ref. jf|. However in contrast to the 
method reported there which can only be applied to classical Hamiltonian without linear 
term in classical momentum, our method can be applied to such a case. 

Moreover, in this specific case where Q is separable Q(q, A; t) = p(q; t)P(X) and P(X) is 
given by Eq. (jSJ), Eq. ( l22|) becomes 

d tP + d qi (g ij (d qj S Q -A j )p)=0. (27) 
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We can then identify an "effective" velocity field given by 

vl = g^(d qj S Q -A 3 ). (28) 

Since Sq(q; t) is just the phase of the Schrodinger wave function, it turns out that the above 
effective velocity field is equal to the "actual" velocity of particle (beable) in pilot-wave 
theory 6|. This allows us to draw a conclusion that our model will reproduce the prediction 
of pilot-wave theory on statistical wave-like interference pattern in slits experiment and 
tunneling over potential barrier However unlike pilot-wave theory, in our model, the 
dynamics is strictly stochastic, the wave function is not physically real-field and the Born's 
statistical interpretation of wave function is valid for all time by construction. 



III. MEASUREMENT WITHOUT WAVE FUNCTION COLLAPSE AND EXTER- 
NAL OBSERVER 

Now let us apply the quantization method developed in the previous Section to a class 
of von Neumann measurement model. To do this, let us consider the dynamics of two inter- 
acting particles, the first particle with coordinate q\ represents the system whose properties 
being measured, and the second particle with coordinate q 2 represents the measuring ap- 
paratus. Let us first discuss how classical dynamics describes the whole system. Let us 
suppose that one wants to measure a quantity A x = A^qijp ). To do this, let us choose the 
following classical measurement-interaction Hamiltonian: 

K = gA 1 (q 1 ,p i )p 2 , (29) 

where g is a coupling constant. Let us further assume that the interaction is impulsive so 
that the individual free Hamiltonians of the particles are ignorable. 

In classical mechanics we first have dA x jdt = {A x , H} = 0. A x is thus conserved during 
the measurement-interaction. The idea is then to correlate the value of A^qxjP ) with the 
classical momentum of the apparatus p 2 while keeping the value of A x {qx,pS) remained un- 
changed. On the other hand, one also has dq 2 /dt = {q 2 ,H_} — dAx, which can be integrated 
to give 

q 2 (T) = q 2 (0) + gA.T, (30) 

where T is time span of the measurement-interaction. The value of A ± prior to the measure- 
ment can thus be inferred from the observation of the initial and final values of q 2 . In other 
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words, the position of the second particle plays the role of the pointer of the measurement 
apparatus. The above model is of course far from realistic. Especially, our model excludes 
the irreversibility of the registration process which can only be done by realistic apparatus 
plus bath with large degree of freedom. See for example Ref. 8] for a realistic model of 
measurement. Here we shall focus on the issue of wave function collapse. 

Now let us consider an ensemble of identically prepared system whose classical Hamil- 
tonian is given by Eq. (129]) . and quantize it as in the previous section. For concreteness, 
without loosing generality, let us consider measurement of angular momentum. To make 
explicit the three dimensional nature of the problem, let us put qi = (xi,yi,Zi). First let us 
consider the measurement of z— part angular momentum of the first particle 

L Z1 = xiP yi - yip xi , (31) 

where p x is the conjugate momentum of x\ and so on. The measurement-interaction classical 
Hamiltonian of Eq. ( 129]) thus reads 

Ki = gL zl p 2 = g(xiP yi - vip Xi )p 2 - (32) 

The Hamilton- Jacobi equation of (EJ) then becomes 

d t S + g{x 1 d yi S - yid Xl S)d q2 S = 0. (33) 

On the other hand, substituting Eq. (132]) into Eq. fl5]), the classical velocity field is given by 

xi(S) = -gyid q2 S, y^S) = gxid q2 S, z^S) = 0, q 2 (S) = g(x 1 d yi S - yid xl S). (34) 

The continuity equation of (jl]) then becomes 

dtp - gyid xi (pd q2 S) + gxid yi (pd g2 S) + gx x d q2 (pd yi S) - gyid q2 (pd xl S) = 0. (35) 

Next, keeping in mind Eq. (134]) . the rule of replacement of Eq. ([7]) thus becomes 

P^Q, d Xl S^d Xl S + ±^, d yi S^d m S+^, d q2 S^d q2 S+^, 

d t S H- d t S + + g\{ Xl d yx d q2 S - yid Xl d q2 S). (36) 

Let us proceed to see how Eq. ( 136]) modifies Eqs. ( 133]) and ( 135]) . Imposing the first four 
equations of (136]) into Eq. (135]) one obtains, after a simple calculation 

d t n - gyxd^Qd^S) + gx^iVld^S) + gx^^ttd^S) 

-gyid q2 (Vtd Xl S) - gX(yxd xl d q2 Q - xid yi d q2 Q) = 0. (37) 



On the other hand, imposing the last four equations of (|36|) into Eq. (|33|) . one has, after an 
arrangement 

d t S + g(x l d yi S - yi d Xl S)d Q2 S - g\ 2 - y^^-) 

+ WFi( dt ^ ~ gyid Xl (ttd q2 S) + gxxdy^Vtd^S) + gxid q2 (Qd Vl S) - gyid q2 (Vtd Xl S) 



2Q 

-gX(y 1 d Xl d Q2 Q - xid yi d q2 VL) ) = 0, (38) 



where we have again defined a real- valued function R = yVt and used the identity of Eq. 
(|T8|) . Substituting Eq. (|37|) into Eq. (1381) . the last term in the bracket vanishes to give 

d t S + g(x 1 d yi S-y 1 d Xl S)d q2 S-gX 2 [x 1 ^^-y 1 ^^) = 0. (39) 



The dynamics of ensemble of trajectories is then determined by pair of coupled Eqs. (J3j 
and ( l39l) which are parameterized by the hidden random variable A. 

Again, assuming the symmetry of Eq. (I2"0"j) . one can see that S(q,X;t) and S(q,—X;t) 
satisfy the same differential equation of (13"9"1) . Hence, assuming that initially S(q, A; 0) = 
S(q, — A; 0), one obtains S(q, A; t) = S(q, — A; t). As in the previous section, this can then be 
used to eliminate the last term on the left hand side of Eq. (|37|) to give 

d t tt - gyi d xl (Qd q2 S) + gx 1 d yi (ttd q2 S) + gx 1 d g2 (nd yi S) - gyi d q2 (Qd Xl S) = 0. (40) 



Further, defining complex- valued wave function as in Eq. ( 1231) . Eqs. ( 1391) and ( 1401) can then 
be rewritten into the following modified Schrodinger equation: 

A 2 

i\\\d t q = -gX 2 {x l d yi - yi d Xl )d q2 ^ = g—L zl p 2 V, (41) 

where L Zl = —ih(xid yi —yid Xl ) and p 2 = —ihd q2 are the quantum mechanical z— angular mo- 
mentum and linear momentum operators pertaining to the wave functions of the first and sec- 
ond particle, respectively, and again we have assumed that the spatiotemporal fluctuations 
of A is ignorable as compared to that of S. Finally, assuming that f2(g, A; t) = p(q; t)P(X) is 
separable and taking the case when -P(A) is given by Eq. ([H]), the above modified Schrodinger 
equation reduces into 

iW t ^ Q = H&q, with H t = gL Zl p 2 , (42) 

where ^Q(q;t) is the Schrodinger wave function defined as in Eq. ( )25|) . This result can be 
extended to the measurement of angular momentum along the x— and y— directions by 
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cyclic permutation of (x,y,z). In this case, L Zl in Eq. (J12j) is replaced by L Xl and L yi , 
the quantum mechanical angular momentum operators along the x— and y— directions, 
respectively. We have thus reproduced the results of canonical quantization as a specific 
case of our hidden variable model. Further, in this case, Eq. fHOj) becomes 

dtp - gyid xl (pd q2 S Q ) + gxid yi (pd q2 S Q ) + gxid q2 (pd yi S Q ) - gyxd^pd^Sq) = 0, (43) 

from which we can extract an effective velocity field which is equal to the actual velocity 
field of particles in pilot-wave theory. This then allows us to follow all the argumentation 
of pilot-wave theory to describe measurement without wave function collapse and external 
observer reproducing the prediction of quantum mechanics. 

As is discussed in Ref. [3j, the above method can also be applied to measurement of 
position and momentum. Let us remark however that unlike measurement of linear and 
angular momentum, the measurement of position is special. Namely, we can show that in 
this case the pair of classical quantities {p, S} and the corresponding quantum quantities 
{p, Sq} satisfy the same Schrodinger equation with quantum Hamiltonian H = gq\p2- Hence, 
there is no quantum correction to the original classical equations. In other words, unlike 
measurement of linear and angular momentum, the measurement of position can reveal 
the value of the position of the particle prior to measurement. This implies important 
consequence. Recall that the results of measurement of linear and angular momentum (or 
measurement of any physical quantities) are inferred from the position of the second particle, 
namely the apparatus pointer (in reality, any model of measurement should be reducible to 
position measurement). Then one might argue that one needs another particle, the third 
particle, as the second apparatus to probe the position of the second particle (the first 
apparatus). Proceeding in this way thus will lead to infinite regression: one will further 
need the forth particle (the third apparatus) to probe the position of the third particle (the 
second apparatus) and so on. In our model, however, since the quantum treatment of the 
position measurement is equivalent to the classical treatment revealing the position of the 
particle prior-measurement, then the second measurement on the position of the second 
particle (the first apparatus) is in principle not necessary. Namely, the results of position 
measurement by the second, third, forth apparatuses and so on are all equal to each other. 
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IV. POSSIBLE CORRECTION TO QUANTUM MECHANICAL PREDICTION 



We have shown in the previous two Sections that the prediction of quantum mechanics 
is reproduced corresponding to a specific distribution of hidden random variable A given by 
Eq. (jSJ). It is then imperative to see the implications if one allows |A| to fluctuate around h 
with small yet finite width. In this case, instead of using Eqs. ( 125]) or (|42|) . one has to work 
with Eqs. ff24|) or f l4Tj) . and regards them as the natural generalization of the Schrodinger 
equation. 

Let us first discuss measurement of angular momentum in ensemble of identically prepared 
system so that the initial wave function of the system (first particle) ip(qi) is given by one of 
the eigenfunction of the angular momentum operator tp(q\) = <fii{qi), L Zl <pi = Itfii, where I is 
the eigenvalue. Let us denote the initial wave function of the apparatus (second particle) by 
V^ofe), assumed to be sufficiently localized. The total initial wave function of the system- 
apparatus is thus given by 



We have thus made an idealization that the initial wave function is independent of A. Recall 
that in this case, according to the standard quantum mechanics, each single measurement 
event will give outcome I with certainty (probability one). This is one of the postulate of 
quantum mechanics in addition to the Schrodinger equation. 

Let us solve Eq. (jUJ) with the initial condition given by Eq. (jUJ). To do this, let us 
assume that after time-span t of measurement-interaction, the wave function can be written 
as ^(q, A;£) = 4>i(qi)f(q2, A; t). Inserting this into Eq. (j4Tj) and keeping in mind that 
L Zl 4>i = l<pu one h as 9 t (f + gl'd q2 (f = 0, with /'(A) = ^-l. It can then be directly integrated 
with the initial condition <^(g 2 , A; 0) = y^ofe) to give <£>(g 2 , A;t) = y?o(<?2 — gl't). One thus 
finally has \l/(g, A; t) = <t>i(qi)tpo{q2 — g\X\lt/K). Hence, in each single measurement event, the 
wave function of the apparatus becomes correlated to the initial state of the system and is 
shifted an amount of gl'(X)t. This means that at the end of each single measurement event, 
the initial position of the second particle (the apparatus pointer) is shifted uniformly and 
randomly as 



Now let us interpret the above formalism in similar way as with classical measurement. 
As discussed in the previous section, in the latter case, after time-span of measurement- 




(44) 



q 2 (t,\) = q 2 (0) + gl'(X)t. 



(45) 
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interaction t, the position of the apparatus-particle is shifted as qz(t) = 92(0) + gL zl t. 
From this, one infers that the result of measurement to be given by L z . Similarly, it is 
natural to interpret Eq. (I45p that the outcome of each single measurement event is given 
by /'(A) = \\\l/h. Hence, instead of obtaining a sharp value / as postulated by the standard 
quantum mechanics, one obtains a random value Z'(A) which depends on the value of the 
hidden variable A. One can also see that when the distribution of A is given by Eq. (jSj) so 
that A = ±h, then the randomness of the outcome of single measurement disappears and one 
regains the prediction of quantum mechanics: l'(±h) = I with probability one. For general 
distribution of A satisfying Eq. 0, we have thus a random correction to the prediction of 
quantum mechanics: even when the initial wave function of the system is given by one of 
the eigenfunction of the angular momentum operator, the result of each single measurement 
will still be random with statistical properties determined by the distribution of A. This 
observation in turn leads to a finite broadening of the spectral line purely induced by the 
ridden random variable. Detail elaborations of this observations is reported somewhere else 



9] 



Next, let us notice that the generalized Schrodinger equation of Eqs. ( |24l) or ( flTi) are still 
linear. Hence, given two solutions ^i(q, X;t) = yQexp(iS/\X\), i = 1,2, one can construct 
new solution through linear superposition = a^i + 6^2, where a and b are complex 
numbers. Calculating the probability density of the position, assuming that f2j is separable 
A; t) = Pi(q; t)P(X), the interference term then takes the following form 

h2(q;t) ~ ^P~ 2 I d\P(\)exp Q| - (46) 

which, for general type of -P(A) = P(—X) and Si(q, A; t), will give different results from the 
quantum mechanics. In particular if -P(A) takes a form of symmetric log-normal function 
then there will be Gaussian suppression of the quantum mechanical interference even in a 
closed system [101 ]. 



V. CONCLUSION AND DISCUSSION 



We have proposed a quantization method based on replacement of c-number by another 
c-number parameterized with hidden random variable. The results of canonical quantization 
is reproduced if the hidden random variable A can only take discrete values ±h with equal 
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probability. Of course A can be a function of other continuous hidden random variables. 
Unlike canonical quantization, the method is free from operator ordering ambiguity and has 
a straightforward interpretation as statistical modification of classical dynamics of ensemble 
of trajectories. Hence, the quantum-classical correspondence is kept physically transparent 
in the quantization processes. The classical limit of the Schrodinger equation is given by the 
classical dynamics of ensemble of trajectories. 

In particular, for all the system considered, we can identity an "effective" velocity field 
which turns out to be equal to the "actual" velocity of the particle in pilot-wave theory. This 
then allows us to conclude that the model reproduces the statistical wave-like interference 
pattern in slits experiment; and further we can borrow all argumentation of pilot- wave theory 
on measurement without wave function collapse and external classical observer. However 
unlike pilot-wave theory, the model is stochastic, the wave function is not physically real 
and the Born's statistics is valid for all time by construction. Moreover, the construction is 
unique given the classical Lagrangian or Hamiltonian. Finally, assuming that |A| fluctuates 
around h with a very small yet finite width, then the model predicts small correction to the 
prediction of quantum mechanics. This might lead to precision test of quantum mechanics 
against our hidden variable model. 

It is then imperative to ask how our model will deal with Bell's no-go theory. Since our 
model reproduces the prediction of quantum mechanics for specific distribution of A, then 
for this case, it must violate Bell inequality which implies that it is non-local in the sense of 
Bell [11 1, or there is no global Kolmogorovian space which covers all the probability spaces 
of the incompatible measurement in EPR-type of experiments 12] , or both. We believe that 
this question can be discussed only if we know the physical origin of the the general rules of 
replacement postulated in Eq. (j7]). To this end, a discussion on the derivation of the rules 
from Hamilton- Jacobi theory with a random constraint is given some where else jl3 |. 
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